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DESIGN STUDY OF SHAFT FACE SEAL WITH SELF-ACTING LIFT AUGMENTATION 

I - SELF-ACTING PAD GEOMETRY 
by John Zuk, Lawrence P. Ludwig, and Robert L. Johnson 
Lewis Research Center 

SUMMARY 

A parametric design study was made on a face seal with a self-acting pad geometry 
having a 6. 5-inch (16. 5-cm) nominal diameter. The load carrying capacity of the self- 
acting pad geometry was calculated for various pad recess depths, film thicknesses, 
angular deformations, and self-acting pad geometries. The self-acting pad geometries 
had a characteristically steep gradient (high gas film stiffness) of lift force against film 
thickness. This high gas film stiffness is advantageous in that the lift force drops off 
rapidly if the seal tends to open; hence, the closing force increases rapidly. If the seal 
tends to close, the opposite effect occurs. A small operating gas film thickness, there- 
fore, is maintained, and the high gas film stiffness forces the seal nosepiece to dynam- 
ically track the runout motion of the seal face. The shallower recess depths produce 
steeper gradients of lift force against film thickness, but for wear considerations, a 
deeper recess depth of 0.001 inch (0.003 cm) was selected. Angular deformation of the 
sealing face did not have a significant effect on the load carrying capacity of the self- 
acting pads. 


INTRODUCTION 

Shaft seal systems in advanced turbine engines will be operated at speeds, temper- 
atures, and pressures higher than those currently used. Conventional face seals pres- 
ently in use in some gas turbine engines are limited to a sliding velocity of about 350 feet 
per second (107 m/sec), a gas temperature of 800° F (700 K), and a. pressure of 125 psi 
(86. 1 N/cm ) (ref. 1). Advanced engines will require seals to operate at 500 feet per 
second (152 m/sec) (ref. 2), and pressures to 500 psi (344 N/cm ) (ref. 3) are antici- 
pated. Face seals for these high speeds and pressures must operate with positive seal- 
ing face separation (no rubbing contact) in order to achieve long life and reliability. Some 



conventional face seals used in current engines apparently operate with positive face 
separation under steady-state conditions, but this positive face Separation is not assured 
by the design procedure. Extension of the speed, pressure, and temperature limitations 
of a conventional face seal will be minimal unless the design ensures the small but vital 
positive sealing face separation. 

A direct approach to obtaining positive sealing face separation is found in the hydro- 
static and self-acting seal designs under investigation in the studies reported in refer- 
ences 4 and 5. The self-acting seal design, which has a low leakage potential (ref. 5), is 
a conventional face seal with a self-acting pad geometry (step-type gas -lubricated bear- 
ings) for lift augmentation (see fig. 1). It should be noted that, in the conventional face 
seal, divergent sealing faces lead to loss of seal opening force; hence, rubbing contact 
can occur (ref. 5). However, as pointed out in reference 5, the step-type gas bearing 
added to a conventional face seal acts to maintain sealing surface separation even though 
adverse sealing face deformation occurs. The feasibility of the self-acting face seal has 
been demonstrated in several NASA sponsored programs (ref. 5). Recent data show 
satisfactory performance to a speed of 400 feet per second (122 m/sec), a pressure of 
250 psi (172 N/cm^), and a sealed gas temperature of 1000° F (810 K). 

The design of the self-acting geometry is vital to seal dynamic performance. The 
self-acting geometry must provide 

(1) Sufficient lift force to separate the surfaces at less than idle speed 

(2) Sufficient gas film stiffness to make the nosepiece dynamically track the axial 

motions of the seal seat face 

A necessary step in optimization of lift augmentation in a face seal is a parametric 
study of the self-acting geometry. Thus, the objectives of the present work are to deter- 
mine for a seal of practical interest 

(1) A practical pad recess depth for the self-acting geometry 

(2) The effect of the number of self-acting pads 

(3) The effect of thermal deformation on load capacity 

The parametric study was made on a face seal with a self-acting lift geometry having 

a 6. 5-inch (16. 5-cm) nominal diameter. Sealed gas pressures of 65, 215, and 315 psia 

2 

(45, 148, and 217 N/cm abs) were considered. The load carrying capacity of the self- 
acting geometry was calculated for various pad recess depths, film thicknesses, angular 
deformations, and for a number of self-acting pads. 


DESCRIPTION OF FACE SEAL WITH SELF-ACTING GEOMETRY 

Figure 1 shows a face-type seal with a self-acting pad geometry that consists of a 
series of shallow recesses arranged circumferentially around the seal under the sealing 


2 



Sump 



dam. It should be noted that the lift geometries are bounded at the inside and outside 
diameters by the sealed pressure P^. This is accomplished by feed slots that connect 
the annular groove directly under the sealing dam with the sealed pressure in the cavity. 
The effect of the self-acting geometry on face seal operation is illustrated in figure 2, 
which shows parallel sealing faces operating without rubbing contact because of a balance 
between the closing force and the opening force, that is, the self-acting lift force plus the 
pressure acting between the sealing faces. If the seal tends to close, the gas bearing 
force increases to prevent rubbing contact; thus, a condition of no rubbing contact can 
prevail except at startup and shutdown. 
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Nosepiece position 




(b) Seal with exaggerated face deformation (coning or dishing). 


Figure 2. - Self-acting face seal. 



RESULTS AND DISCUSSION 


Design Points 

For the purpose of analysis, three design points were selected as representative of 
operating conditions in the advanced gas turbine (see table I). 


Model of Self-Acting Lift Geometry 

A single self-acting pad (shrouded step-type), which served as the mathematical 
model, is shown in figure 3. The radial width b of all pads was 0. 20 inch (0. 507 cm), 
and the length c was determined by the number of pads arranged circumferentially under 
the sealing dam. Because of the large radius-to-pad width ratio r/b, the curvature was 
neglected in the mathematical model. Thus, the pad boundaries conform to a Cartesian 
coordinate system. The feed-groove width f was held constant for all pads. The recess 
length-to-land length ratio R/L was varied between 0.4 and 1. 8. Also, the effect on 
load capacity of the ratio of recess depth to film thickness A/h m was determined over 
the range of (A + h m )/h m = 0. 52 to 6. 0. (All symbols are defined in appendix A. ) Cal- 
culations of load capacity were made for the parallel film face and for the case with a 
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Figure 3. - Model of single self-acting pad (shrouded Rayleigh step type). 


5 



2-milliradian angular deformation. This angular deformation was chosen as a practical 
case that might occur as the result of thermal gradients. Figure 4 shows the actual 
dimensions of one of the pad designs. 



Figure 4. - Details of sealing dam and self-acting pad geometry. All dimensions 
are in inches (cm). 


The computer program listing in reference 5, which is valid for a shrouded lift pad, 
is used for this study. This computer program has been modified (as shown in appendix B) 
so that it can be used on the Lewis Research Center computer, and the input has been 
modified for NAMELIST to facilitate parametric design studies. The program solves the 
complete two-dimensional compressible flow Reynolds lubrication equation written in the 
P 2 form 


where 


J_/h^ 9P 2 \ 9 /h 3 3P 2 \_ A 9(PH) + A a (PH) 

ax \2 ax/ ay\2 ay/ x ax y ay 
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The formulation of this analysis is given in reference 6. For details of the solution of 
the foregoing equation, see appendix B in reference 7. A discussion of the use of the 
program for a nonparallel film is given in appendix A of reference 5. The following as- 
sumptions or restrictions apply in the present analysis: 

(1) The fluid is Newtonian and viscous. A laminar flow regime is assumed. It should 
be noted that it is possible for the seal assembly to be operated in the transition or slip 
flow regime for clearances less than 0.0001 inch (0.0003 cm), and the analysis is not 
valid in this slip regime . 

(2) The bulk modulus, 3X + 2/u = 0, is Stokes idealization and is valid since local 
shocks are not present and the flow is analyzed for the continuum flow regimes. 

(3) The body forces are negligible . 

(4) The flow is laminar for a maximum film thickness of 0.001 inch (0.003 cm). The 
maximum Reynolds number is 658, which is considerably less than the minimum transi- 
tion rotational Reynolds number of about 1900. 

(5) The modified or reduced pressure flow Reynolds number is much less than 1. 

For example, at design point 3 in table I, the modified Reynolds number is in the range of 
0. 023. Thus, the Reynolds equation should be valid for the model of the self-acting lift 
pad. 


TABLE I. - DESIGN POINTS 




Design Point 



1 

2 

3 

Velocity, V, ft/sec (m/sec) 

200 ( 61 ) 

500 (153) 

450 (137) 

Sealed gas pressure, P^, 

65 (45) 

215 (148) 

315 (217) 

psia (N/cm 2 abs) 




Sealed gas temperature, 

100 (311) 

800 (700) 

1300 (977) 

Tj, °F (K) 




Sump pressure, P2, psia 

15 (10.3) 

15 (10. 3) 

15 (10. 3) 

(N/cm 2 abs) 


4. 87x10"® (3. 36X10" 10 ) 

6 . 00x10“® (4.14X10" 10 ) 

9 

Viscosity, \i , (lb)(sec)/in. 

2. 75x10"® (19xl0~ 10 ) 

((N)(sec)/ cm 2 ) 




Operating condition 

Descent 

Takeoff and climb 

Cruise 


The compressibility number A y was calculated for each of the three design points 
(see table I) at various film thicknesses (see table H) by using the following relations: 
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Design point 1: 
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Effect of Number of Self-Acting Lift Pads 


The load capacity of the self-acting pads as a function of the number of pads is shown 
in figure 5. As the number of pads increases, the circumferential length of each pad 


Oj re. 



Figure 5. - Load capacity of shrouded self-acting lift pads 
as function of number of pads. Sealed pressure, 

315 psia (217 N/cm 2 abs) ; sealed gas temperature, 1300° F 
(977 K) ; sliding velocity, 450 feet per second (137 m/sec). 
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necessarily decreases in order to fit the available region under the sealing dam. The cal- 
culations were made at design point 3, which has a pressure differential of 300 psi (206 
N/cra ) and a sealed gas temperature of 1300° F (977 K). The sliding velocity was 
450 feet per second (137 m/sec) .. At a film thickness ratio of (A + h m )/h m = 2, which is 
a practical ratio for a large operating film thickness, an optimum is indicated at 28 pads. 
However, at a film thickness ratio of (A + h m )/h m = 6, which is a practical ratio for a 
small operating thickness, no optimum is indicated and the load capacity increases as the 
number of pads decreases. Since the seal is expected to operate at various film thick- 
nesses, a compromise in the number of pads must be sought, and the number selected 
was 20. 


Effect of Recess Length -to-Land Length Ratio 

The load capacity as a function of recess length-to-land length ratio R/L for sym- 
metric and asymmetric designs is shown in figure 6. The calculation was made for de- 
sign point 3 of table I with the film thickness ratio of (A + h m )/h m = 6. A symmetric de- 
sign has a slightly greater load capacity than an asymmetric design, and the curve trend 
suggests that a recessed land length ratio R/L of 1. 8 or greater will provide near opti- 
mum capacity. Because of wear considerations, it is desirable to have the land length as 
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Figure 6. - Load capacity of shrouded self-acting lift pads 
as function of recess length-to-land length ratio. Sealed 
pressure, 315 psia (217 N/cm 2 abs) ; sealed gas tempera- 
ture, 1300° F (977 K); sliding velocity, 450 feet per 
second (137 m/sec); number of pads, 20; film thickness 
ratio, 6; compressibility number in circumferential 
direction, 1.78. 
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long as practical; therefore, a recess-to-land length ratio of 1.4 was selected as a com- 
promise between wear and load capacity considerations. 



Effect of Recess Depth-to-Film Thickness Ratio 

The load capacity as a function of film thickness ratio (A + h m )/ h m is shown in fig- 
ure 7. The optimum film thickness ratio is near (A + h m )/h m = 2. However, under dy- 
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Figure 7. -Load capacity of shrouded self-acting lift pads as function 
of film thickness ratio. Sealed pressure, 315 psia (217 N/cm 2 abs); 
sealed gas temperature, 1300° F (977 K); sliding velocity, 450 feet 
per second (137 m/sec); number of pads, 20. 


namic operation, the film thickness is expected to vary because of the inherent runout of 
the seal seat face. Thus, in operation, a range of film thickness ratios is expected. 
Because of leakage considerations, the film thickness of 0.0004 inch (0.0010 cm) is a 
practical average operating condition; and for this case, the recess depth would be 
0.0004 inch (0.0010 cm) to achieve the optimum value of (A + h m )/h m = 2. However, in 
order to provide increased wear capability, a recess depth of 0.001 inch (0.00254 cm) 
was selected. Thus, for a film thickness of 0.0004 inch (0.0010 cm), the film thickness 
ratio is (A + h m )/h m = 3.5. 
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Load Capacity as Function of Film Thickness 

The load capacity for the three design points listed in table I is shown in figure 8. 


Design Sealed pressure^ Temperature, 
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Figure 8. - Load capacity of shrouded self-acting lift pads as 
function of film thickness. Number of pads, 20; recess length 
to-land length ratio, 1.4. 


The calculations are for 20 pads and a recess depth of 0. 001 inch (0. 003 cm). The curves 
reveal a high film stiffness that is advantageous for seal operations; that is, at low film 
thicknesses, a high lift force is produced to prevent rubbing contact and at high film 
thicknesses, which would have high leakage, only a low lift force is produced. Thus, if 
the seal opens, the pad lift force drops off sharply, and the closing force acts to return 
the nose to the equilibrium condition. 


11 



Effect of Face Angular Deformation 


The effect of face angular deformation on load capacity is shown in figure 9. A de- 
formation of 2 milliradians was selected as a practical value that might occur as a result 
of axial thermal gradients in the seal faces. This angular deformation results in a loss 
in load capacity. For example, figure 9 shows a 27-pound (12-kg) load capacity for par- 
allel faces and a 21-pound (9. 5-kg) load capacity for faces with a mean film thickness of 
0.0004 inch (0.0010 cm) and a 2-milliradian deformation. This loss in load capacity is 
not excessive and indicates that the self-acting geometry can accommodate some face 



Film thickness, h m , in. 



Figure 9. - Load capacity of self-acting geometry as function of 
angular deformation. Sealed pressure, 315 psia (217 N/cm 2 
abs>; sealed gas temperature, 1300° F (977 K); sliding velocity 
400 feet per second (122 m/sec); number of pads, 20 ; recess 
length-to-land length ratio, 1.4; recess depths, 0.0010 inch 
(0. 0025 cm) and 0. 0004 inch (0. 0010 cm). 
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deformation. (See ref. 5 for discussion of techniques in minimizing deformations in 
seals. ) 


Effect of Recess Depth on Film Stiffness 

The 0. 001-inch (0. 003-cm) recess depth was selected as a practical value because a 
compromise between wear and load capacity was necessary. A comparison of the 0. 001- 
inch (0. 003-cm) depth with the 0. 0004-inch (0. 0010-cm) depth demonstrates the penalty 
paid for the compromise. Figure 10 shows load capacity as a function of film thickness 
for each of the two recess depths. The 0.001-inch (0. 003-cm) recess depth has a greater 
load capacity except at very small film thicknesses, that is, at 0. 0002 inch (0. 0005 cm) 
and below. The important point is that the shallow recess has a higher film stiffness 



Film thickness, h m , cm 

Figure 10. - Effect of recess depth on load capacity of self-acting 
lift pads. Sealed pressure, 315 psia (217 N/cm? abs) ; sealed 
gas temperature, 1300° F {977 K); sliding velocity, 450 feet 
per second (137 m/sec); number of pads, 20; recess length-to- 
land length ratio, 1.4. 





(curve of steeper load against film height), and this higher film stiffness is desirable in 
seal operation; that is, if the seal tends to open, the self-acting force drops off rapidly 
and, hence, the closing force increases rapidly. If the seal tends to close, the opposite 
effect occurs. 


SUMMARY OF RESULTS 

A parametric study was made on a face seal with a self-acting lift geometry having 

a 6. 5-inch (16. 5-cm) nominal diameter. Sealed gas pressures of 65, 215, and 315 psia 

o 

(45, 148, and 217 N/cm abs) were considered. The load carrying capacity of the self- 
acting geometry was calculated for various recess depths, film thicknesses, angular 
deformations, and numbers of self-acting pads. A self-acting geometry was placed on a 
high-pressure side of the seal dam inside diameter. The analysis yielded the following 
results: 

1. The characteristic steep gradient of lift force against gap height of the self-acting 
geometry was responsible for maintaining the small sealing gap height necessary for low 
leakage. That is, if the seal tends to open, the lift force drops off rapidly; hence, the 
closing force increases rapidly. If the seal tends to close, the opposite effect occurs. 
Thus, the self-acting geometry provides a high gas film stiffness, which is necessary if 
the nosepiece is to track dynamically the face runout of the seal seat. 

2. The gradient of self-acting lift force against film thickness was affected by the 
recess pad depth; the shallower recesses produced a steeper gradient. 

3. For faces with angular deformation (e.g. , due to thermal gradients), the self- 
acting lift force was less than that for parallel faces. Calculations for practical values 
of angular deformation indicated that the self-acting geometry produces a useful load. 

4. The self-acting pad geometry selected was a compromise between wear and load 
capacity considerations; that is, the recess depth selected was deeper than the optimum, 
and the recess length-to-land length ratio was smaller than the optimum. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 7, 1970, 

126-15. 
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APPENDIX A 


SYMBOLS 

2 2 

A pad area, in. ; cm 

b self-acting pad radial width, in. ; cm 

c self-acting pad circumferential length, in. ; cm 

f feed-groove circumferential length, in. ; cm 

H recess depth plus film thickness, in. ; cm 

h film thickness, in. ; cm 

L land circumferential length, in. ; cm 

o 

P static pressure, psi; N/cm 

R recess circumferential length, in. ; cm 

r radius, in. ; cm 

T temperature, °F; K 

U moving seal seat radial surface speed, ft/sec; m/sec 

V moving seal seat circumferential surface speed, ft/sec; m/sec 

W load capacity, lb; kg 

X radial coordinate direction 

Y circumferential coordinate direction 

a relative inclination angle of surfaces, mrad 

A recess or pocket depth, in. ; cm 

2 

A„ compressibility number in radial direction, 6/iUc/P 0 h 
x » in 

2 

Ay compressibility number in circumferential direction, 6juVb/P a h m 

X second viscosity coefficient or coefficient of bulk viscosity 

2 2 

ju absolute or dynamic viscosity, (lb)(sec)/in. ; (N)(sec)/cm 

Subscripts: 

a ambient conditions 

m mean 

1 sealed conditions 

2 sump conditions 


APPENDIX B 


COMPUTER PROGRAM FOR SHROUDED SELF-ACTING LIFT PAD 

Input Instructions 

The computer program is composed of the main program and five subroutines. The 
main program reads the input, performs the major portion of the calculations, calls on 
the subroutines, and writes the output. The five subroutines are titled SUM, HFUN, 
HXFUN, HYFUN, and MATIN V. 

SUM is the Simpson’s rule integration scheme for odd and even mesh points. 

HFUN describes the shape of film thickness. 

HXFUN describes the derivative of film thickness in X-direction. 

HYFUN describes the derivative of film thickness in Y-direction. 

MATINV inverts the necessary matrix required by the columnwise influence coeffi- 
cient method. 

The program is coded in FORTRAN IV, version 13, and is run on an IBM 7044-7094 
direct couple system with 32K core storage. 

Each set of input data requires three major input cards. In addition, the first set 
requires the seal geometry cards. All input statements are in the NAMELIST form, a 
feature of FORTRAN IV, version 13. The following is a list of the data cards for a sample 
problem: 

(1) Pad Number Control Card 

Item Number (3 items) 

(a) NPAD-Control number specifying number of pads to be used 

1 - Single pad 

2 - Twin pad, extra card (lb) is needed 

3 - Twin pad, NPD = TRUE, extra card (la) is needed 

- Triple pad, NPD = FALSE, extra cards (la and lb) are needed 

(b) INP-Control for input to be continued or terminated 

0 - Input set to be continued 

1 - Input set to be terminated 

(c) NDPP-Control for printout W/A (Pg - Pj) 

0 - W/A (Pg - Pj) will be printed out besides W/A Pj 

1 - Only W/A Pj printed out as the load parameter 
(la) Second Pad Geometry Card 

Item Number (6 items) 

1 - LA, number of mesh points to first step in X-direction for second pad 

2 - JA, number of mesh points to first step in Y-direction for second pad 
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3 - IAA, number of mesh points to second step in X-direction for second pad 

4 - JAA, number of mesh points to second step in Y-direction for second pad 

5 - STEDA, ratio of step depth to film thickness for second pad 

6 - NPD = .TRUE, for twin pad 

= .FALSE, for triple pad, additional card (lb) is required 
(lb) Second or Third Pad Geometry Card 
This card is required when NPAD = 2 or NPAD = 3 and NPD = FALSE , and it 
describes the second or third pad geometry, respectively. 

Item Number (5 items) 

1 - IB, number of mesh points to first step in X-direction for second or 

third pad 

2 - JB, number of mesh points to first step in Y-direction for second or 

third pad 

3 - IBB, number of mesh points to second step in X-direction for second or 

third pad 

4 - JBB, number of mesh points to second step in Y-direction for second or 

third pad 

5 - STEDB, ratio of step depth to film thickness for second or third pad 
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Figure 11. - Sample mesh -point layout. (All nodal points without an assigned value 
should be set equal to zero.) 


(2) Geometry, General Control, and Film Geometry Card 
Item Number (29 items, see fig. 11) 

1 - M, number of mesh points in X-direction (13 points max) 

2 - N, number of mesh points in Y-direction (30 points max) 

3 - PLAMX = 6/iUb/P a h^, bearing number in X-direction 

9 

4 - PLAMY = 6juVb/P a h m , bearing number in Y-direction 
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5 - YOX, circumferential length (Y-direction)/radial width (X-direction) for 

a single pad 

6 - IH, number of points to first step in X-direction for first pad 

7 - JH, number of points to first step in Y-direction for first pad 

8 - IHH, number of points to second step in X-direction for first pad 

9 - JHH, number of points to second step in Y-direction for first pad 

10 - STEDE, ratio of step depth to film thickness for first pad 

11 - NDIG, number of digits required for accuracy, usually for single pad set 

12 - PFIX (1), dimensionless pressure, P^/Pj = 1 

13 - PFIX (2), dimensionless pressure, P 2 /Pj 

14 - PFIX (3), dimensionless pressure, Pg/P 4 

15 - NCASE, identify sets of input data 

16 - LKOUNT, number of iterations allowed, usually set LKOUNT = 25 

17 - IFLO = 2, for inlet flow in X-direction 

18 - JFLO = 2, to N-l, for flow in Y-direction 

19 - IFLOE = M-l, for outlet flow in X-direction 

20 - COE (1) = Cj/L 

21 - COE (2) = C 2 /L 

22 - COE (3) = 1 

23 - COE (4) - C 4 /C 3 

24 - COE (5) = C 5 /C 3 

25 - COE (6) = C 6 /Cg 

26 - COE (7) = C 7 /Cg 

27 - COE (8) = Cg/Cg 

28 - COE (9) = Cg/Cg 

29 - COE (10) = C 10 /Cg 

The film thickness can be expressed in the third-degree polynomial as 
h = Cg + C 4 (X-C 1 ) + C 5 (Y-C 2 ) + CgCX-Cj) 2 + C 7 (Y-C 2 ) 2 

+ c 8 (x-c 1 )(y-c 2 ) +C 9 (x-C 1 ) 3 +c 10 (Y-c 2 ) 3 


For parallel film study, COE (3) = 1.0, and the rest of the coefficients are equal to zero. 
As for the detail for pivoting and tilting of the seal pad, reference 5 should be consulted. 
(3) Final Control Card 

Item Number (4 items) 

1 - QUEP1 = T, print out P 2 after each iteration 
F, omit P 2 
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2 - PP0UT1 = T, print out final P after convergence 

F , omit final P^ 

3 - POUT1 = T, print out final pressure distribution 

F, omit final pressure distribution 

4 - NUKUEl = T, read complete set of seal geometry cards, usually used 

for first set 
(4) Seal Geometry Cards 

This set of data cards contains code numbers that describe each of the MXN nodal points. 
Each nodal point is designated as KUE in the computer program language. Since a two- 
dimensional array is stored columnwise in the computing machine, the MXN code num- 
bers must be entered in columnwise order. 

KUE = 0, regular point or corner of depressed area 

KUE = 1, 2, or 3, corresponding to given pressure PFIX (1, 2, or 3) 

KUE = 4, vertical line of step in X-direction 
KUE = 5, horizontal line of step in Y-direction 
KUE = 6, TOP JOINT, leading edge in X-direction 
KUE = 7, BOTTOM JOINT, trailing edge in X-direction 
KUE = 8, LEFT JOINT, leading edge in Y-direction 
KUE = 9, RIGHT JOINT, trailing edge in Y-direction 
Note that the joint condition implies that the boundary specified satisfies the cyclic bound- 
ary requirements. At present, the computer program is set up to handle correctly the 
cyclic boundary condition in the Y-direction. For this boundary condition, the mesh 
points in the Y-direction are increased from N to N +1, and N - 1 divisions to N divisions 
for the total length of the seal pad. Furthermore, the input of N for cyclic boundary 
conditions should be equal to N, not N + 1, because the computer program increases the 
mesh points N to N + 1 automatically. The calculations for flow rate in the Y-direction 
for cyclic boundary conditions should set JFLO = N. 

Comments on Use of Shrouded Rayleigh-Step Computer Program 

The following remarks should be useful to a potential user of this program. Some of 
the limitations are stated. 

(1) The number of mesh points is limited to 30 points in the Y-direction and 13 points 
maximum in the X-direction because of the 32K core storage capacity of the Lewis 
Research Center computer. Thus, a difficulty arises in checking to see if there is solu- 
tion convergence of a particular geometric configuration and operating condition. 

(2) The convergence of the numerical analysis is limited to moderate compressibility 
numbers. This difficulty may be overcome by incorporating the Newton -Raphson proced- 
ure. It is believed that convergence at higher compressibility numbers will be possible. 
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(31 An error occurs in the analysis in the step region where large gradients occur. 
This may explain the behavior shown in figure 6, where an optimum pocket length-to-land 
length ratio is not found. 


Program Listing 


SIEFTC TRIPAD LIST, DECK 

R I NC 10 REVISED BY DR.CHAS.NG 7,11,67 FOR MULTIPLE DEPRESSED FEED ANO 
C. CORRECTED THE INTEGRATION FCR CYCLIC BOUNDARY CONDITION IN Y DIRECTION 

C PROGRAM TO SOLVE STEP COMPRESS. BEAR. PROB. WITH FIXED 

C BOUNDARIES, LINES OF SYMMETRY, JOINTS IN ANY DIRECTION. 

C EQUATIONS ARE WRITTEN FOR PARALLEL FACES ONLY 

C M , N ARE THE MESH POINTS IN THE X AND Y DIRECTIONS 

C KUE=0 REGULAR POINT OR CORNER OF DEPRESSED AREA OR LINE OF SUMM 

C KUE = 1 * 2 , 3 KNOWN PRESSURE= PFIX(1,2,0R 3) 

C KUE = 4 VERTICAL LINE OF STEP 

C KUE=S HORIZONTAL LINE OF STEP 

C KUE=6 TOP JOINT 

C KUE = 7 BOTTOM JOINT 

C KUE=B LEFT JOINT 

C KUE = 9 RIGHT JOINT 

C PROBLEM IS SOLVED COLUMNWISE. M (FIRST INDEX) 

C SHOULD BE SMALLER THAN N(SECCNC INDEX! 

C X IN I DIRECTION (VERT. DOWN) 

C Y IN J DIRECTION (HOR. LEFT TO RIGHT) 

C PLAKX,PLAMY=X,Y, COMPONENTS OF PLAM 

C ( IH, JH ) , ( IHH, JHH) ARE CORNERS OF STOP BOUNDARY 

C STEDF=STEP DEPTH. WHERE NO STEP H=1 

C ( IB, JB ) , ( IBB, JB8 ) ARE CORNERS FOR SECOND PAD OR POCKET 

C STEDB, STEP DEPTH FOR 2 PAD OR POCKET 

C ( I A, JA) , ( IAA, JAA) ARE CORNERS FOR SECOND OR THIRD PAD ANO POCKET 

C STEDA,STEP DEPTH FOR 2 OR 3 PAD AND POCKET 

C ND I G= NO OF DIGITS WANTED REPEATED TO TRUNCATE SOLUTION 

C LKOUNT IS THE MAXIMUM ALLOWABLE NUMBER OF ITERATIONS 

C I F L 0= I COORDINATE OF THE LINE ACROSS WHICH Y-FLOW IS COMPUTED. 

C JFLO= J COORDINATE OF THE LINE ACROSS WHICH X-FLOW IS CGMPUTED. 

C COE= CLEARANCE COEFF I C I EN T S . SEE HFUN, HXFUN, HYFUN. 

C QREP=. TRUE. PUT OUT P2 AFTEREACH ITERATION 

C PPQUT= . TRUE . OUTPUT OF P2 AFTER CONVERGENCE 

C POUT=. TRUE. WANTED OUTPUT CF P. AFTER CONVERGENCE 

C NE WKUE = . TRUE . IF NEW KUE ARRAY IS READ IN 

C TWO PAD SET NPD= TRUE AND THREE PAD SET NPD = F AL SE 

C NPAD=NO.OF PADS, SINGLE PAD=1,TW0 PAD=2, THREE PAD=3 

C NDPP = 0 W/A/DP PRINTOUT,. USE CNLY WHEN P2 IS GREATER THAN PI 

C NDPP^l USE ONLY FOR REGULAR GAS BEARING STUDIES 

COMMON E , R , D 

DIMENSION P F I X ( 3) ,KUE( 13,3C),QFIX(3),H(13,30), 
1Q(13,30),FF(13,30)»F(13),A(13,13),R{13),C(61),E(13,13,.30), 

1QSMA( U, 13) ,G( 13,30) ,R( 13, 13,30 ) ,S( 13,30) ,PP (61 ) ,PX(61 ) , 

IPY(61 ) ,D< 13, 13, 30) ,CQQ( 61 > ,COE (10) 

1 * X X ( 13) , Y Y ( 6 1 ) ,PFX( 13,30) »PFY( 13,30),HFX(13,30) ,HFY( 13,30) ,F5( 13, 

130 ) 

DIMENSION LOGCON ( 4 ) * LOGKN 1(4) 

DATA L0GTRU/1HT/ 
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EQUIVALENCE { LQGCON { 1} , QREP ) , ( LOGCON ( 2 ) ,PPCUT) ,! LOGCON ( 3 ), POUT ) , 

A ( LOGCON ( 4 ) , NEWKUE ) , ( LCGKN1 ( 1 ) , QREP l ) , (LOGKN1 (2) ,PPOUT 1 ) » 

B { LOGKN 1(3) , POUT 1 ) , < LOGKN1 (4) .NUKUEl ) 

LOGICAL JOINT, QREP , PPOU T , PCU T , NEWKUE , NPO 
LOGICAL LCGCON 

NAMELIST /CARD1/NPAD, INP.NDPP 

NAMELIST /CARDIA/IA, JA, IAA, JAA,STEDA,NPD 

NAMELIST /CARDIB/IB, JB, I BB , JBB , STEDB 

NAMELIST /C ARD2/ M, N, PL AMX, PLAMY, YOX, IH, JH, IHH, JHH, STEOE , NO I G , PF IX, 
A NCASE ,LKOUNT, I FLO, JFLO, IFLCE,COE 
NAMELIST /C ARD3/ QREP 1,P POUT I , POU T l , NUKUE 1 
NAMELIST /CARD4/KUE 
1 FORMAT ( IX70I I ) 

3 FORMAT! 25H MATRIX IS SINGULAR AT J= I3,16H,CASE A B ANOONED . / IH 1 ) 

A FORMAT 1 / 10! IX, FI 1.7) ) 

5 FORMAT! //18H CASE CONVERGES TO I3,14H DIGITS AFTER I3.11H I TER AT I 
IONS ) 

6 FORMAT ( // 23H FINAL RESULTS FOR CASE I5//13H FORCE/AREA =E14.7, 

1 64HC0EF OF CENTER OF PRESSURE IN PERCENTAGE OF SIDE 
2D I MENS I UN S = (E14«7»1H» E14.7.2H).) 

7 FORMAT (46HOFLOW PER UNIT LENGTH IN X AT ENT. AND EXll =( 

1 1PE 14. 7 , 2H, 1PE14.7 , IH) 24HF LOW PER U. L. IN Y = ( 1PE14.7, IH )/ ) 

f< FORMAT ( 29H1F INAL PRESSURE DISTRIBUTION. //) 

9 FORMAT ( 25H1F INAL P**2 DISTRIBUTION. /) 

90 2 FORMAT! lH0,3HIA=,l5,3HJA=,I5,4HIAA=,I5,4HJAA=,l5,6HSTEDA=,E15k8) 

90 3 F0RMAT(IHC,3HIB=,I5,3HJB=,I5,4HI8B=,I5,4HJBB=,I5,6FSTEDB=,E15.8) 
905 FORMAT ( // 23H FINAL RESULTS FOR CASE I5//13H FORCE/AREA =F14.7,7HW/ 
1A/0P=»E14.7/65H- COEF OF CENTER OF PRESSURE IN PERCENTAGE OF SICE 
2DIMENSI0NS = (E14.7,1H, E14.7,2HI.) 

12 READ ( 5 , C AR D 1 ) 

GO TO ( LO ,707, 701 ) , NPAD 

701 READ! 5,CARD1A> 

IF ( NPD ) GO TO 10 
707 READ! 5 ,CARD1B) 

10 READ ( 5 , C AR02 ) 

READ ( 5 , C ARD3 ) 

WRITE (6,1002) QREPL ,PP0UT1 ,P0UT1 , NUKUEl 

1002 FORMAT (32H1 QREP PPOUT POUT NEWKUE/6XA l ,6XA l , 8XA1 ,6XA l ) 

DO 1004 M M = l ,4 

IF ( LOGKN 1 ( MM ) .FQ. LOG TRU ) GO TO 1006 
LOGCON ( MM ) = .FALSE. 

GO TO 1004 

1006 LOGCON ( MM ) = .TRUE. 

1004 CONTINUE 

WRITE (6,1003) M , N , PLAMX, PLAMY, YOX, IH, JH, IHH, JHH, STEDE, 

2 NO I G , PFIX, NCASE, LKOUNT, IFLO, JFLO, I FLOE , COE 

1003 FORMAT ( IH0,6X, IHM, 5X,3HN ,21HPLAMX PLAMY YOX , 7 X , 2H I H , 6 X , 2H JH 

2 5 X , 3H I HH 5X » 3HJHH/1H I 7 , I 6 , 3F8 . 3 , 4 1 8 / 80H0 STEDE NO I G PFIX! 1 ) 

3 P F l X ( 2 ) PF I X ( 3 ) NCASE LKOUNT IFLO JFLO IFL0E/1H F7.3,.I 
45, F10. 3, 2F8. 3, 17, 218, 19, 17/6 HO COE = 10F8.3) 

GO TO 1702,703,704) , NP AD 

704 WRITE (6,902) I A , J A , I A A , J A A , S T E D A 
IF ( NPD ) GO TO 702 

703 WRITE (6,903) I B , J B , I BB , JBB , ST EDB 

702 IF! .NOT. NEWKUE )G0 TO 35 
20 READ l 5,CARD4) 

DO 30 1=1, M 

WRITE! 6, l ) (KUE! I ,J) ,J=1,N) 

DO 30 J= 1 ,N 


30 KUE ( I , J )=KUE ( I , J ) + l 
35 K0UNT=0 
NN=N— 1 
MM=M— 1 

DO 40 K = 1 , 3 

40 QF I X < K ) = PF I X ( K ) *PF I X ( K ) 

DX=1./FLGAT (MM) 

IF(KUE(2,1).EQ.9.AND.PLAMY.GT.PLAMX) GO TO 44 
DY=YOX/FLCAT (NN) 

GO TO 43 

44 DY=YOX/FLCAT { N ) 

43 00 41 1 = 1 » M 

4 L XX ( I ) = F L 0 A T ( I-l)*DX 

00 4? J=i , N 
4 2 YY ( J ) = F LO AT ( J - 1 )*DY 
S A = 1 • / D X 
SB= 1 • / D Y 
S A A = S A* S A 
S8B=$B*SB 
SC=SAA+SBB 
SD=-2 . *SC 
SE=PL AMX / ( 2 . *DX ) 

SEE=1./(2.*DX) 

SG=PLAMY/ (2.*DY) 

SGG=1./(2.*DY) 

SH = 2.*SA 
SI=2.*SS 
S J = ?.* SAA 
SK=2.*S88 
DO 50 1=1, M 
DO 50 J=1,N 

5 0 H ( I,J)=HFl'N(XX(I ) f Y Y ( J ), f COE ) 

GO TO (51,52,53) , NP AD 
53 00 61 I = I A , I A A 
DO 61 J = J A , J A A 

61 H( I , J )=H( I , J ) + S TEDA 
IF ( NPD ) GO TO 51 

52 DO 6? i = I B , IBB 
DO 62 J = J R , JBB 

62 HU , J )=H< I , J ) + S TEDB 
51 DO 60 1= IH, IHH 

DO 60 J = J H , J HH 
60 H ( I , J)=H< I , J ) + S T EDE 
DO 70 1=1, M 
DO 70 J = l , N 
0 ( I , J ) = 1 

HFX ( I , J ) =3.*HXFUN( XX ( I ) ,YY ( J ) ,COE ) /H ( I , J ) 

HFYU , J ) =3.*HYFUN( XX ( I ) , YY ( J ) ,COE ) / H ( 1 , J ) 
ZZZ=-i.O/(H( I,J)*H( I, J) ) 

F5( I ,J ) = ( 2./ 3. )*ZZZ«*(HFX ( I , J ) *PLAMX+HFY ( I , J )*PLAMY ) 
PF X ( I , J)=PLAMX*Z?Z 
70 PF Y ( I , J ) =PLAMY*ZZZ 
JO INT=. FALSE. 

8 <J DO 90 1=2, MM 

IF (KUE (1,1) .EQ. 9* OR. KUE ( I , 1 ) .EG. 10 ) JOINT=. TRUE. 

90 CONTINUE 
100 DO 130 I = l , M 
G( I, 1 )=0. 

DO 110 K = 1 , M 
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110 E(I »Kf L ) = 0 • 

IF ( .NUT. JOINT) GO TO 130 

00 1?0 K— 1 1 M 
0( I * K , 1 ) = 0. 

1 F ( I . E(J . K ) 0 ( I,K,L)=L.O 
120 CONTINUE 

130 CONTINUE 

DO 370 J=1,N 
00 310 I-1,M 

F F ( I ,.J ) = 1.0/ SORT ( ABStQt I,J) )) 

Ft I ) =0 • 

KU“ KUE ( I t J ) 

GO T0( 140,210,210,210,230,250,270,270, 140, 1 AO > ,KU 
140 SGGG = SGG *t FF t I , J) *PFY ( I , J ) +HFY ( I , J > ) 

C ( I ) = SBB + SGGG 
B( I) = S BB- SGGG 
00 150 K -= 1 , M 
A( I , K ) =0 • 0 

IFtK.FQ. [ )A( I,K)=SD+F5( ItJ)*FFII,J) 

150 CONTINUE 

IF(JOINT) GO TO 180 
IF ( J.NE. 1 ) GO TO 160 
8 ( I ) =0 . 

Ct I )=SK 

160 IF (J.NF.N) GO TO 170 
C ( I ) =0 . 

B ( I ) =SK 

170 IF < I .EQ. 1 > GO TO 190 
IF(I.EQ.M) GO TO 200 
180 SR = SEE*(FF( I ,J )*PFXt I , J)+HFX( I ,J) ) 

At I , 1 + 1 )=5 SAA+SR 
At I» I - 1 ) =SAA-SR 
GO TO 310 
190 At I , H-1) = SJ 
GO TO 310 
200 At l , 1-1 ) = SJ 
GO TO 310 

210 KK U = KUE t I , J )-l 
Bt I ) =0 . 

Ct I >=0. 

Ft I ) = QF I X < KK U ) 

DO 220 K = 1 , M 
At I f K)=0. 

IF ( I . EQ . K ) A ( I ,K ) = l .0 
220 CONTINUE 
GO TO 310 

230 HP LUS = H ( I ,J) 

HMINUS^Ht I , J ) 

GO TO (231,232,233) tNPAD 
233 IF (J.FU.JA) HMINUS=HMINUS-STECA 
IF (J.EQ.JAA) HPLUS=HPLUS-STEDA 
IF t NP 0 ) GO TO 231 

232 IF (J.EQ.JB) HM l NUS-HM I NUS- S TE CB 
IF (J.FU.JBB) HPLUS=HPLUS-STEDR 

231 IF (J.EQ.JH) HMINUS-HMINUS-STEOE 
IF (J.EQ.JHH) HPLUS=HPLUS-STEDE 
HHH=HPLUS**3 

HH=HMINUS**3 

S8=.5*SB 




B ( I ) = —SB * HH 
C ( I ) = — S8 * HHH 
DO 240 K= 1 , M 
240 A ( I » K ) = 0. 

A ( I » I ) = S8 * ( HH+HHH ) + PL AM Y* ( HPLUS— HMINUS ) +FF ( I, J) 
GO TO 810 

250 HPLUS=H ( I , J ) 

HM I NUS = H ( I , J ) 

IF (IIH .EQ. IB) .OR. (IH .EG. I A > ) GO TO 251 
GO TO (25 1,252,253) ,NPAD 
253 IF(I.EO.IA) HMINUS=HMINUS-STEDA 
IF (I.EO.IAA) HPLUS=HPLUS-STEOA 
IF (NPD) GO TO 251 
252 IF(I.EQ.IB) HM INUS=HMINUS-STEDB 
IF (I.EO.IBB) HPLUS=HPLUS-STEDB 

251 IF( I.EQ.IH) HMINUS=HMINUS-STEDE 
IF (I.EQ.IHH) HPLUS=HPLUS— STEOE 
HH = HM I NUS **3 

HHH=HPLUS**3 
B ( I ) =0 . 

C( I ) =0 . 

D9 = 2 • *DX 
DO 260 K= 1 » M 
260 A I I * K ) =0. 0 

A ( I , I ) = ( hH+HHH) /D9-PLAMX#{ HM INUS-HPLUS ) *FF ( I , J ) 
All, 1+1 )=-HHH/D9 
A ( I , I-l )=-HH/D9 
GO TO 310 

270 SGGG=SGG* ( FF (I,J)*PFY(I,J)+HFY(I,J)) 

B ( I ) = SBB- SGGG 
C( I ) =SB8+SGGG 
DO 280 K= l , M 
A ( I , K ) =0 . 0 

IF(K.EQ.I) A ( I ,K )=SD +F 5 ( I , J ) *F F ( I , J ) 

280 CONTINUE 

SR = SEE*(FF( I , J ) * PFX ( I , J ) +HF X ( I,J) ) 

I F ( KU. EQ.8) GO TO 290 
I F ( KU.EQ.7) GO TO 300 
GO TO 310 
290 A ( 1,1) =S4A+SR 
A ( I , I-l ) =SAA-SR 
GO TO 310 
300 A ( I ,M)=SAA-SR 
A ( I, I+1)=SAA+SR 
GO TO 310 
310 CONTINUE 

DO 320 1=1, M 

DO 320 K = 1 » M 

320 QSMA ( I , K ) =A ( I , K ) + B ( I ) *E ( I , K , J ) 

CALL MATINV(QSMA,M,BB,0,DET,ID) 

GO TO (340, 330), ID 
340 DO 360 I = l , M 
G ( I , J + 1 ) = C . 

DO 360 K= 1 » M 

G( I , J+l ) =G( I , J + l) + Q S M A ( I,K) + (F(K) — B(K)*G(K,J) ) 

E( I ,K, J+l )=-GSMA ( I ,K)*C (K) 

IF ( .NOT. JOINT ) GO TO 360 
DUM=0.0 

DO 350 KK=1,M 
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350 DUM-DUM-QSMA ( I,KK)*B(KK)*D<KK,K,J) 

0 ( I , K , J + 1 ) = D UM 
360 CONTINUE 
370 CONTINUE 
DM A = 0 . 0 

IF(JOINT) GO TO 410 
DO 330 1 = i,M 

DM A = AM AX 1 (DMA, ABS(Q( I »N)-G( I tN+1) ) ) 

380 Q ( I , N ) =G ( I ,N+1 ) 

DO 400 J J =2 » N 
J = N + 2- J J 
DO 400 I = 1 , M 
DUM=0. 0 
DO 390 K = 1 , M 

390 DUM = DUM + E ( I , K, J ) *Q( K, J ) 

DUM=DUM+G( I , J) 

DM A = AM AX 1 (OMA, ADS(DUM-G( I , J-L ) ) ) 

400 Q ( I , J-1)=DUM 
GO TP 560 
410 DO 420 1= 1 , M 

DO 420 K=1,M 
QSM A ( I t K ) = — D ( I ,K,N+1 ) 

IF ( I.EO.K ) QSM A ( I , K ) = QS M A { I ,K) + 1.0 
420 CONTINUE 

CALL MATINV(QSMA,M,RB,0,DET,ID) 

GO TP (430,330) , ID 
330 WR I T E ( 6 t 3 ) J 
GO TP 12 

430 DO 460 1 = 1, M 

DU = 0.0 
DO 450 K = 1 , M 
DUM= 0.0 
DO 440 KK = 1 , M 

440 DUM=DUM+QSMA( I ,KK)*E(KK, K , N+l ) 

R( I ,K,N)=OUM 

450 DU = DU + t)SMA( I ,K)*G(K,N+1 > 

460 S( I ,N) =DU 

DO 490 J J = 2 , N 
J = N+2-JJ 
DO 490 1=1, M 

DU = 0 • 0 

DO 480 K = 1 , M 
DU M = 0.0 
DO 470 K K = 1 , M 

470 DUM= D( l ,KK,J)*R(KK,K,N)+E( I,KK,J)*R(KK,K, J)+0UM 
R ( I , K , J - 1 ) = DUM 

480 DU=DU+D( I,K,J)*S(K,N)+E(I,K,J)*S(K,J) 

490 Si I , J-l ) = CU + G ( I , J) 

DM A = 0 . 0 

DO 500 1 = 1, M 

DO 500 K= 1 , M 

QS M A ( I , K ) =-R ( I , K , 1 ) 

IF( I.EQ.K)QSMA( I ,K)= QSMA(I,K)+1.0 
500 CONTINUE 

CALL MAT INV(QSMA,M,BR,C,DET,ID) 

GO TD(51 0,330) , ID 
510 DO 530 1 = 1, M 
DU = 0 . 0 

DO 520 K = 1 , M 
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520 DU=DU+QSMA( I ,K)*S(K, I ) 

DM A = AM AX i ( DMA, ABS(DU-Q< 1,1))) 

530 0(1, 1 )=DU 

DO 550 J= 2 , N 
DO 550 I - 1 , M 
DU = 0 • 0 

DO 540 K=1,M 

540 DU = DU+R< I ,K, J)*Q(K, 1 ) 

DU = DU+ S ( I , J ) 

DMA = AMAX1 (DMA, ABS(DU-Q( I , J ) ) ) 

550 Q( I , J )=DU 

560 IF(QREP) WRITE(6,4)Q 
KOUNT = KOUNT+ I 

I F ( KDUNT • GE . LKQUNT ) GO TO 561 
I F ( DM A .GT. 10.0**FL0AT(-NUIG) )G0 TO ICO 

561 WRI ] F ( 6 , 5 )NDIG,K PUN T 
IF(PP0UT)feRITE<6,9) 

DO 575 1 = 1, M 

IF { PP OUT ) WR ITE ( 6,4) (Q( I , J ) , J=1 ,N ) 

DO 570 J = 1,N 

570 Q( I,J)=SQRT(ABS(0(I, J) ) ) 

575 CONTINUE 
IF(P0UT)WRITE(6,8) 

IF(KUE(2, 1).EQ.9.AND.PLAMY.GT.PLAMX) GO TO 571 
DO 576 1=1, M 

IF ( POUT )WRITE(6,4)(Q <I,J)tJ=l,N) 

576 CONTINUE 

GO rn 573 

571 NE =N + 1 

DO 572 1 = 1, M 
H ( I , NE ) = H ( I , 1 ) 

Q( I , NE ) = Q ! 1,1) 

IF (POUT) kR ITE (6,4) ( Q ( I , J ) , J= 1 , NE ) 

572 CONTINUE 

573 IF (KUE (2, 1 ) •EQ.9.AND.PLAMY.GT.PLAMX ) GO TO 7C0 
DO 590 1=1, M 

X= DX*FL 0 A T ( I-l ) 

DO 5R0 J = 1 , N 
5 8 C QOQ ( J ) =U ( I, J )-l.O 

PP(I)= SUM ( QCQ , N , DY ) 

590 PX ( I )= PP ( I )*X 
DO 610 J = 1 , N 
Y = FLOAT! J-1)*DY 
DO 600 1 = 1 ,M 
600 OQQ ( I )=Q( I , J >-1.0 
610 PY(J) = SUM(CQQ,M,r>X)*Y 
F P = SUM ( P P , M , DX ) 

F X = SUM(PX,M,DX) 

XF •= F X / F P 
F Y = SUM ( P Y , N , D Y ) 

GO TP 745 
700 DO 710 1 = 1, M 

x=dx*float ( i-i ) 

DO 720 J = 1 , N E 
720 OQQ! J )=Q( I , J ) - l • 

PP ( I ) - SUM (QQQ, NE,DY ) 

710 PX! I )=PP( I )*X 
DO 730 J= 1 » NE 
Y=FLPAT ! J-l ) *DY 
DO 740 1=1, M 



7*0 ^00 < I )=w( It J )-L. 

730 PY(J)=SUM(C0Q,M,OX)*Y 
FP = SU*MPP,M,DX) 

F X = SUM ( PX ,M , UX ) 

XF = FX/FP 
F Y = SUM ( P Y , NE , DY ) 

745 YF=FY/FP/Y0X 
FP=FP/ YOX 

IF (NDPP) 8 80, 8 BO ,881 

880 FDP=FP/(PFIX(2)~PFIX(1) ) 

WRITF (6,905) NC ASE , FP , FDP , XF , YF 
GO TO 882 

881 WR I T F ( 6 , 6 ) NC ASE , F P , XF , YF 

882 DO 620 1=1, M 

62 0 PP (I)=Q< I « JFLO)*H( I 9 JF LO ) * { -PL A MY + H ( I , JFLO ) **2 
1 * ( Q ( I ,JFLC+1 )-Q( ItJFLO— 1) ) / ( 2 - * D Y ) ) 

FLOY = SUM( P P , M , DX ) 

NFL □= 1 
IFLT=IFLO 

625 IF ( KUF ( 2 , 1 ) . EQ.9. AND • P L AM Y . G T • PL AMX ) GO TO 750 
DO 630 J = 1 , N 

630 PP ( J ) = Q{ I FLT, J >*H( IFLT, J ) * ( -PL A M X +H ( I FL T , J ) **2 
1 * ( 0 ( IFLT+1,J) — Q{ IFLT-lrJ) ) / ( 2 • * D X ) ) 

FLOT=SUM( PP,N,DY) /YOX 
GO TO 770 

750 DO 760 J = 1 , iN E 

76 0 PP< J)=U(IFLT, J )*H{ IFLT* J)*(-PLAMX+H( IFLT, J >**2 
1 *( Q ( IFLT+lfJ) — CM IFLT— 1».J) ) / ( 2 • *DX ) ) 

FL 0 T= SUM ( P P , N E , D Y ) / YO X 

770 GO TO ( 635,636 ) , NF L 0 

635 NF L 0 = 2 
FLOX=FLOT 
IFLT=IFLOE 
GO TO 625 

636 CONTINUE 

WR I T F ( 6 , 7 ) FLOT , FLOX , FLOY 
I F ( INP ) 12, 12,800 

800 STOP 
END 


$ I BF TC SUM LIST, DECK 

FUNCTION SUM(P,M,DX) 
DIMENSION P(70) 

I F ( (M/2)*2-M) 80;81,81 

80 K= 2 

KK=M-1 
KKK=2 
SUM=0. 0 

10 DO 20 I =K , KK , KKK 
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20 SUM=SUM+P(I) 

GO TO (30,40,50) ,K 
30 SUM=SUM*DX/3.0 
RETURN 
40 K= 3 

45 SUM = 2.0*SL'M 
GO TO 10 
5 0 K=L 
KK = M 
KKK=M-1 
GO TO 45 

81 SUM = P ( 1 ) — PCM) 

SUM = SUM + . 2 5*P ( 1 ) + 3 . *P ( 2)-.25*P(3) 
DO 60 J=3,M,2 

SUM = SUM +4.*P ( J ) +2.*P ( J + l ) 

60 CONTINUE 

SUM = SUM *DX/3. 

RETURN 

END 


tIBFTC U FUN LIST, DECK 

FUNCTION HFUN <XX,YY,COE) 

C FUNCTION TO EVALUATE CLEARANCE H FROM X,Y AND COEFFICIENTS COE 

DIMENSION COE ( 10 ) 

X= XX-COE ( 1 ) 

Y= YY-COE ( 2 ) 

HFUN=CQE (3)+C0E(4)*X+C0E(5)*Y+C0E(6)*X*X+ 

ICOE ( 7)*Y*Y + C0E ( 8 )*X*Y + COE ( 9 ) * X * X *X + 

2C0E ( 10 ) *Y*Y*Y 
RETURN 
END 


$IBFTC HXFUN LIST, DECK 

FUNCTION HXFIJN( XX, YY ,COE ) 

C FUNCTION TO EVALUATE X-OERIVATIVE OF H FROM X,Y ANC COEFFICIENTS C 

DIMENSION CUE ( 10) 

X = X X-C OF ( 1 ) 

Y= YY-COE ( 2 ) 

HX FUN-CQE (4 ) +2.*C0E (6)*X+CCE(8)*Y+ 

13.+C0E ( 9 ) *X*X 
RETURN 
END 
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$ I BF TC HYFUN LIST, DECK 

FUNCTION HYFUN(XX,YY,CCE) 

C FUNCTION TO EVALUATE Y-DERIVATIVE OF H FROM X,Y AND COEFFICIENTS C 

DIMENSION COE ! 10) 

X= XX — COE ( 1) 

Y= Y Y — COE ( 2) 

HYF UN = COE ( 5 ) +2 * *COE ( 7) * Y + COE { 8 ) *X+ 3 * *COE ( 10>*Y*Y 

RETURN 

END 


SI8FTC MAT INV LIST, DECK 

MATRIX INVERSION WITH ACCOMPANY I NG SOLUTION OF LINEAR EQUATIONS 
NOVEMBER 1962 S GOOD DAVID TAYLOR MODEL BASIN AM M AT 1 
SUBROUTINE MATIN V (A, N1,B, Ml, DETERM, ID) 

C GENERAL FORM OF DIMENSIGN STATEMENT 

DIMENSION A ( 13, 13 ) ,8 ( 1 3, 1 ) , INDEX! 13 ,3) 

EQUIVALENCE ( I ROW, JROW) , ( I COL CM , JC OL UM ) , ( AMAX, T , SWAP) 

M= M l 
N = N l 

10 DE TE RM = 0 . 0 
IS DO 20 J=1 , N 
20 INDEX! J , 3 ) = 0 
30 DO 650 1 = 1, N 

C SEARCH FOR PIVOT ELEMENT 

40 AM A X = 0 . 0 
45 DO 105 J = 1 , N 

I F ( INDEX ( J , 3 )- 1 ) 60, 1C5, 6C 

60 DO 100 K = 1 , N 

IF(INDEX(K,3)-1) 80, 100, 715 

80 IF ( AMAX -ARS (A(J,K))) 85, 100, 100 

86 I RO W= J 
90 I C OL UM = K 

AMAX = ABS ( A( J,K) ) 

LOO CONTINUE 
105 CONTINUE 

INDEX! I CO LUM , 3 ) = INDEX! IC0LUM,3) +1 

260 INDEX! I , 1 ) = IR0W 
270 INDEX! I , 2 ) = I COLU M 

C INTERCHANGE ROWS TO PUT PIVCT ELEMENT ON DIAGONAL 

130 IF ! IROW- ICOLUM ) 140, 310, 140 

140 DETER M=- DETERM 
150 DO 200 L = 1 , N 
160 S W A P = A ( I R 0 W , L ) 

170 A! I ROW , L )=A( ICOLUM, L) 

200 A ( ICOLUM , L ) = SW AP 

I F ( M ) 310, 310, 210 

210 DO 250 L = 1 , M 
220 SWAP = B( I R C W , L ) 

230 B! IRPW.L ) =B< ICOLUM, L) 

250 R ( ICOLUM, L)=SWAP 


c 


DIVIDE PIVOT ROW BV PIVOT ELEMENT 
310 PIVOT =A{ ICOLUM, ICOLUM) 
DETERM=DETERM*PI VO T 
330 A( ICOLUM, ICOLUM) = I.O 
34 0 DO 3 SO L = 1 » N 

350 4(ICnLUM,L)=A(IC0LUM,L)/PIVCT 
35 S I F ( M ) 380, 380, 360 
360 DO 370 L = 1 , M 

370 R( ICOLUM, L)=B( ICOLUM , L ) /PI VCT 
C REDUCE NON-PIVOT ROWS 

380 DO 5S0 L 1 = 1 , N 
390 IF (Ll-ICOLUM ) 400, 550, 400 
400 T=A(L1, ICOLUM) 

420 A ( L 1 , ICOLUM ) =0.0 
430 DO 450 L = l , N 

450 A(LI,L)=A(L1,L)-A(IC0LUM,L)*T 
455 I F ( M ) 550, 550, 460 

460 DO 500 L=L,M 

500 B(LI,L)=BIL1,L)-B( ICOLUM, L)*T 
550 CONTINUE 

C INTERCHANGE COLUMNS 

600 DO 710 1=1, N 
610 L = N+ l- I 

620 IF ( I NDE X I L , 1)*— INDEX (L,2) ) 630, 710, 63C 
630 JROW= I IMDEX ( L , 1 ) 

640 JCOLUM= INDEX ( L ,2 ) 

650 DO 705 K = 1 , N 
660 SW AP= A { K , JRQW ) 

670 AIK, JROW) = AIK, JCOLUM) 

700 AIK, JCOLUM)=SWAP 
705 CONTINUE 
710 CONTINUE 

DO 730 K = 1 , N 

IF I INDEX ( K, 3 ) -1) 7 15,720,715 

72 0 CONTINUE 
730 CONTINUE 
I D — 1 

GO TO 740 
715 ID =2 
740 RETURN 
END 


Sample Problem 

2 

An example of the use of the computer program is given for the 315 psia (217 N/cm ) 
design point for both a parallel film and a film with a 2-milliradian tilt. The geometry is 
shown in figures 4 and 11. The mean film thickness to be investigated h m is 0.0003 inch 
(0.0008 cm), and the groove depth A is 0.0001 inch (0.0003 cm). 

From table II for h m = 0. 0003 inch (0. 0008 cm), 


A = 1.715 = PLAMY PLAMX = 0 = A 

y x 


30 



I 


TABLE H. - SUMMARY OF COMPRESSIBILITY NUMBER FOR EACH 
FILM THICKNESS AT EACH DESIGN POINT 


Film thickness , 
m 

Design point 

1 

2 

3 

in. 

cm 

Sealed pressure, 



65 psia (45 N/cm^ abs) 

215 psia (148 N/cm 2 abs) 

315 psia (217 N/cm 2 abs) 




Compressibility number, A 

y 

0.0001 

0.0003 

15.2 

20.4 

15.4 

.0002 

.0005 

3.80 

5.10 

3.86 

. 0003 

.0008 

1.69 

2.27 

1.715 | 

.0004 

.0010 

.950 

1.28 

.965 

.0005 

.0013 

.650 

.816 

.617 

.0006 

.0015 

.423 

. 567 

.428 

. 0007 

.0018 

.311 

.417 

. 315 

.0008 

.0020 

.238 

. 319 

.241 | 

.0009 

.0023 

. 188 

.252 

. 191 

.0010 

.0025 

. 152 

.204 

. 154 


YOX = - = 7 . § = 3. 500 

b 0.25 


STEDE = 


A + h, 


m 


‘m 


3. 3 


PFIX(l) 



1 


PFIX(2) 



1 


For the parallel film case , 



= 0 


= 0 


I 


I 


C 3 = Coe(3) = 1 


c 4 = 0 


For the 2-milliradian tilt case, 


C , = Coe(l) = — = = 0.5 

1 b b 


c 2 = o 


Cg - 1 


C = __ (0- 002)(0. 2 5) _ ^ 00fj 

4 li 0.0003 

m 


The data for this sample problem are shown in tables III and IV. The execution time 
for this problem was about 0. 21 minute on the Lewis Research Center computer. (Note, 
however, that the execution time increases with since the total number of iterations 
increases.) The sample output is as follows: 

For the parallel film case, 


w = force = 0> 03813818 
area 


or 


load = _W =WPbc = 2.517 l b 
pad pad a pad 

load = _W_ = w p be = h 144 k S 
pad pad a pad 


The total load for 20 pads is 50. 34 pounds (22. 90 kg). 
For the 2-milliradian tilt case, 


w = force = 0.03056596 
area 


W = W P„bc 

d 


2.0175 lb 
pad 


32 



The total load is 40. 35 pounds (21. 85 kg). 

A check was made to see if the solution was the correct converged solution by chang- 
ing the mesh size, as shown with the following results: 


M 

N 

Force/area 

7 

14 

0.03863340 

7 

27 

.03758219 

13 

14 

. 03973810 

13 

27 

.03813818 


Since the results agree satisfactorily, the solution is considered converged. 





TABLE IV. - DATA FOR SAMPLE PROBLEM - CASE II 


[Relative inclination angle of surfaces, 0.002 rad.] 


§ FORTRAN STATEMENT 

2 5 U 5 6 T 8 9 10 ii i? 13 li. 15 16 i» if >9 20 ?' ?? '*» 7;. 25 29 V 3? H <*. 5 s ) JftFv S8 IQ 1.0 l, ■ t. ? i ^ ^ uU u \ » ' 

,$ CARPI NP A.p Ul, IN P - 1 f N D P'P = 1 > 


U.9 SO 5i 52 53 Si. SS S6 5? 58 S') 60>«l 62 63 66 65 66.67 60 69 70 71 72 .73 76 75 76 77 78 79 0 


$ t C A RD2 M = 13r,_N = 27 , PLA MX = 0. , PL A MY = 1. 7 15^ Y 0X - 5 .5 00, IH = 5 , JH=1 , I H H=ll, JHH=19|, 


4 t S^TED|E = 3 ( . 4 30> NTI G=j5 , P FI X=1 . 0, , 1^. 0 , ( 0 . 0,NC AS E|= 2,, 


IP\Lj0E==12,, C$E = O t ,5,0.0., 1.0,, ,1 . 6.67, 


$ CARD 3 QREP'l = f F r , P P0 U T 1 = 1 F r , P0UT = 1 T ' 


CARD4 KUE = 13*1, 


LK.jfUN.T =,50 , IF L, 0^ 2, 


T™ - 


1., 0, . 5, , , 7 *0 . 5jJD., 0., b.,0., 1 , 1,, 0., .5., ,7 * 0 . 5, . 0. , 1 ■ 1 . 0 . 5 , ..7,» , 0„ S , , 0,, , 1 „ 


(F o ur copies | of the foregoi n g < dat a car d ar e r e gui red|. ) < 


1, 0, 5 , 7* 0, 5 , 0, 1 , 1 , 0 , 0, 1; 4^ 0, (), 1 


1;, 11* 0., 1, 1, 11*0 11* 0 , , 1, 1., 11*0^1, 1;, 11* O f jl f 1,1 1*0 ,1,15*1 


>2 3 6 5 | 6 | 7 8 9 10 II 12 1 13 U. 15 16 17 ie|l9 20 21 22 23 26 125 26 27 28 29 3o| 31 32 33 36 35 36|37 30 39 60 61 62^3 66 65 66 67 6B|fc9 50 51 52 53 Sfc|s5 56 57 58 59 6o|6l 62 6366 65 66|&7 68 69 70 71 72|?3 76 75 76 77 70 79 I 
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Sample Problem 

Printout 

Case I: Parallel Film 



HKtP r'PLiUi I’OUI .WLrfPtUt. 

f h i r 



.1 U KL.li-U HLAI-IY YuX, 

i-J X 1 0. 1. / IS 5. 50 U 

iri Jri 1HH 

3 1 1 L 

JHH 

19 

3 1 L Ut null; HI- 1 A ( 1 J FhiA(2) PF1A13) 

J.JUu i l.UUu l.udO U. 

filLA5fc LKUUiMT IFLO 

l 50 2 

JFLU IFLUE 

Zb 1 Z 


Cut = U. U. l.duU -u.Oud -u.UOU -U.OdO -u.UuO -0.000 -U.000 -0.000 

UluUJuOujJU LuulUULOUlUOUU 1 
15555 5j3->:>55 

lOOUUUUUUUUU LUUUUu*lULLdUdUl 
luUddUUu Ud 0 U lUdl/UJ iUouUUUU 1 
1 JLULUuudUddLUUUUU'.UUdUOUUi 
1 uuu uuUU LudJLJU uJj^u^uJddu 1 
IOUOUUUOLUULLUg LLU4UGuUOUU1 
luUUUoOdduJj LuUUUjiJOdddddl 
1 LLd uv/JOiyOaU LOUUCJsULUdJUJl 
lj;i3^^a^33jbji^3DdduUUUuUi 
iOLdGUUUUdUuLdUuUUdUUUUUuO i 
iiiiiUilm 111 ill iiliiiiil 


L/Otr Cu.'tVLrloL^ lu J UiUl Tb rth TlK 4 iftKAllUNb 
F iivAL HiMUJiit. 0 1 5 111 ILmJ 1 lUU. 


1. UOUUU uu 
1 . Uuuuuuu 

1 • lUuUbdU 

l.dJdduud 
l.ddUUudd 
1. . OU UU u OU 

1 . OOUuUOu 
1 . UU JUUU J 

i . uuuuudd 

1. UUuOuOU 
1. dddJUJU 
L. UUOdUUU 

1 . UUUJOOO 
l . uOOJUOO 
1 . 0O0J J JO 

i.OOOJOOO 
l.OOOJ JUO 
1 .UOOOdUU 

i.OOOOoOO 

l.OOOOdOO 

1 . oocoaoo 

1.0000000 

1.0000000 

1.0000000 
1. 0000000 

1.0000000 
l. 0000030 

1 . JudddUd 
i •u^obi' >tL 
1- (J 190 J1 / 

L ,\JJ iLtL lv f 
l , \j*. 1 1 0*11 
l.Ul 59* /> 

1 . JdlUOtU 
1 . UJ iZ Ldd 
1 . UUWd/L 

l.dd/lUod 
1. Jiool/ / 
1 . dOo 744 7 

1 . OU 9u Ud 5 
1 . 0900 1 95 
1 . 00 95 93 J 

l.Ul21dU0 
1. U44GUdO 
1.QU22224 

1.0190333 

1.0906100 

1.0000000 

1.01 75609 
1-0511383 

1.02 04360 
1. 0402565 

1.0234120 

1.0284507 

1 . L LUUUUu 

l.ujcuddi 

l.djVuJlo 

1 .UU4/1 /'( 

1.05 o'j'r-v y 

1 . o<; c 5 buj 

1 . UJ9493.' 
1 . JubdddJ 
l. kJlo/Noi 

1 . U 1 9 J 590 

1. o7^uUUo 
1. Ulz9593 

1.0193 l/d 
l .0799 301 
i . UU dJ 00 1 

1.02990 79 
1. 0077022 
1.00 92050 

1.0296967 

1.0950697 

1.00C0000 

1.03 50 594 
1. 1044478 

1.0406700 

1.0783780 

1.0465037 

1.0540475 

1 . LddddOu 
l.u^Vldl 
i.ubloa-to 

A • 00 4 7 H O 0 

i ■ U J *5Uou 
1 .U J CS55u 

l . OU 95 55 a 
1 . UdOudmJ 
1 . \JZ 041 5*1 4 

A . U 144433 
i. 0/ 3040 J 
1. 0 1 l9d 

i . U 199399 
1 . Jb J9910 
1 . 01 1 7 d 90 

1.02436U J 
l.Udb2334 
1.U0 6J 969 

1.0250317 
1.0969 961 
1.0000000 

1.0352774 
1. 1051019 

1.04 09220 
1.1142256 

1.0467909 

1.0742341 

1 . Ll/tlduud 
1 . d5bifi ( 
1 . Uoi't't yv 

l.ddi/od 
1 . U 5 yjujj 
L.dt‘i < tyi i 

A . Uu *5 V4 1 
i . Uu UZ 93 J 
i . u 51 a 5u / 

1. O 1 95 109 
1. 0/35599 
La U£ z 1 y>6 

1 . 0 195 zu3 
1 . JdJ/d *d 
1 .0143923 

1 . 029o 09 6 
1.0b do 13 5 
1.0073916 

1.029964U 
L. 0960693 
1.0000000 

1.0354333 
1. 1055875 

1. 0411023 
1. 1148294 

1.0469964 

1.0844594 

1 . ULuJudU 
i.UbJ^yyo 
l.oobd, Lo 

1 . Ld 4 / c*. o 
a . o? v /2 u7 

nUffllNdV 

1 . UU90£. 5 L 
1 . dddtd f 5 

1. 01 434*3 
1. 0 73 55 J el 

1. JX9UO jli 

1 . U1 95 7d3 
A . Jo 09 97 1 
1 .0 1 oJ 1 75 

1.0297353 
1 .UB009O 1 
1 .0002332 

1.0300935 
1.09 /1269 
1.0000000 

1.0355269 
1. 1058832 

1.0412106 

1.1151813 

1.0471198 

1.0893640 

1 . COdiludu 

1. U5 J / 1 

1 . U 6o 5 S 

A.UU4/0 li 
1 . U 5 'll 7 e. 5 

i.UJLlNC J 

I. dd^udJti 
1 . UU u3 t. 5 d 
1 . U 5 09 50 3 

1. 0 i 43 3 

1 . U /5u 1 J / 
l. u2 34 9t)9 

1 . 0 1 9j 937 
l.UalJuyj 
1 . Oio3u92 

1.02973 73 
1 . Oo 07 23 2 
1 .0005 175 

1.0300700 
1.07 72 126 
1.0000U00 

1.0355582 
1. 1059823 

1.0412467 
1. 11 52971 

1.0471609 

1.0908284 

1 .LOUJUud 
l.Obj^ddc 
i • LoG Ji. Lo 

i.duN/oLd 
l.UJbau/ 
i.utoV'J vo 

1 . ou y bX 3 i 
i. du 090/3 
i.dJbiNtt 

1. U 1 4 3 995 

1 . U / J 33U Z 
1. UXNdOjJ 

1 . Ul93 7o3 
1 .UoU9992 
1 . U 1 Go A 73 

1.0297359 
1. 000096 1 
1.0002332 

1.0300935 
1.07 71269 
1.0000000 

1« 0355270 
1. 1058 632 

1.0412106 

1.1151813 

1.0471198 

1.0893640 

1 . L LUUUUO 

l.dbjli 1 7 
Udoit*.^ 

1 . dd H / L-Jo 

1 • u) IJO’J J 
l.Ut'iibLU 

i. . uu 93 99 1 
I . Uu O/Nll 1 
i.ddldbod 

1. ol 93 105 
1. U 75 5 599 
L. UZ£ I 930 

1 . 0 1 7 3 2 0 3 
1 .Ooo7o7d 
1. 0143923 

1.0296697 
1.0dd6l56 
1 .00 /391b 

1.0259690 
1.09oB694 
1. OUUUOOO 

1.0354333 
1. 1055875 

1.0411023 

1.1148294 

1.0469964 

1.0844595 

1 • L ud J u J u 
J.dbXv 1 Jl 
l.Lbiod'i < 

l . Ud 4 /*, oU 
1 . 0 J IJuuu 
i . u JGV5 J / 

I . 0d /3d3o 

J. . Uu OUUOU 
1 • U*1 Oil 59 5 

1 . 0 A 99933 
L. U / 5093 5 
i. Ol iZ 199 

A . U 1 79 59 9 
1 .00099 1U 
1 .01 1 7 0 9 1 

1.02 95 0 0 3 
1. Ob 02 339 
1. 00 60 96 9 

1.0298317 
1.09 6446 L 
1 . 0000000 

1.0352774 
1. 1051020 

1. 04 0922 0 
l. 1142257 

1.0467909 

1.0742342 

1 . ddddduu 
1. U5236GH 
1.05/0 5 i9 

i.ud't/Wi 
1.1150*44 4 
1 . uZ db 5t>4 

1 . ud 99 95 Z 
1 • UU5 0 UO U 
1 . UA « /4 04 

l.dl^j 390 
1.0 7/dJdO 
A. 0 1 2 9 399 

1.0193 1 7 d 
1 . U /77 30 l 
1 .dddduilX 

1.U299U79 
1. 00 7/022 
1 .00 92 d5U 

1.0296467 

1.U95B647 

l.UOOOOOO 

1.0350594 
1. 1044478 

l. 04 06 700 
1.0783 781 

1.0465037 
L. 0540475 

1 . UuOudUO 
1.0^054.53 
l.dl'Yod'*/ 

L.ddXX /V / 
a . u^: 9 / u*t j 
i . oi 59<r ou 

1 • ud9 009 0 

1 . dJdX loo 
A . UU9 / d7 1 

1. 00 71 009 
1 . 05 00 1 / 7 

i.ddd/'ltd 

l . U09oO do 
1.0900 173 
l .0043433 

1.0121001 
1.09 96 00 1 
I. 0022225 

1.0148333 

1.0486101 

1.0000000 

1.0175610 

1.0511383 

1.0204360 

1.0402566 

1. 0234121 
1.0284507 

1 . uUuuuuo 

1 . UdudddU 
1 . luddddu 

1 . UU LU UUU 
1 . UL UU UU 0 
1 . Uu Lu UuU 

A. udddddd 
1 . UJOuOOU 
1 . ouuuooO 

A. uOUuUUO 
1. UUdJdJU 
1 . JdduddU 

1 . uouJuoo 

1 . OOOJ J JO 
1 . OuUJ 000 

l.OUUOOOO 
i.OOOJOOO 
1. OU0J0OU 

000 
c. c c 

C C. O 

coo 
0 C C 
coo 
0 C O 

1.0000000 
1. 0000000 

1.0000000 

1.0000000 

1.0000000 

1.0000000 

F 1 IN A L K A. juL f o 

l-u* La it 

l 








K» KU/AKL4 = 

J. JJl JdldL" 

JUuir ur LcMlk Uf PKt5 

30 K t in PfcKLlNI AGt JF $1 UE 

UiMENSlUN S 

= ( 0.500000IE 00. 0. 

547521 4E 00) 

FL L>n PtK UI, i 1 

L LHd I H id 

A A I L IN 1 . A'N 0 

C Ai 1 - i-z . 

3d 7 V 1 7 9c. -0 1 « 

2.5d79l6lt-01)FLUH PER U 

. L. IN Y = 

( -1.7639388E 00) 


♦dl» fcA 1 l k.i I iUHaj 
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Case II: Two-Milliradian Tilt 



QREP PPOUT POUT NEHKUE 
F F T I 

H N PLAMX PLAMY YCX, IH JH IHH JHH 

13 27 0. 1.715 3.500 3 1 11 19 

STEDE NDIG PFIXdJ PFIX(2) PFIX<3> NCASE LKOUNT LFLO JFLO IFLOE 
3.300 3 1.000 1.000 0. 2 50 2 26 12 

COE - 0.500 0. 1.000 1.667 -0.000 -0.000 -C.000 -0.000 -O.COO -0.000 

liiuiLiiiiiiii in iiiuiin 

loooooooooooooacooooooooooi 

15555555555 5555 5 5 50000C0001 
100000000000000000400000001 
100000000000000000400000001 
10000000000000QC00400000001 
lOOOOOOOOCOOOOOCOOAOOOOOOOl 
100000000000000C00400000001 
IOOOOOOOCQCOOOOC004000C0001 
100000000000000C00400000001 
155555555555555555000000001 
10000000000000CC00000000001 

iiiiiiiimiiiiiiiiiiiuui 


CASE CONVERGES TO 3 DIGITS AFTER 4 ITERATIONS 


FINAL PRESSURE distribution. 


1.0000000 

1. 0000000 

1.0000000 

1.0000000 
1 .oocoooo 
1.0000000 

l.OCOOOCO 

l.OCOOOCO 

1.0000000 

1.0000000 

1.0000000 

l.OCCOOOO 

l.OCOCCOO 

l.OOOOCOO 

l.OOOCCOO 

ucccoaoo 

uooooooo 

lkCCCOCOO 

1 • CCQCOCO 

ucocooco 
1 * cooco CO 

l.OOOOOCO 

ucooooco 

l.coaoooo 

l.OOOOCOO 

l.ooaoooa 

l.OCCOOOO 

1.0000000 

1.0302417 

1.0608665 

I .0017834 
1.0345379 
1.0504143 

1.00421C8 

1.0391909 

1.04011C2 

1-0068890 

1.0442585 

1.0308482 

1.0097142 

1.0498192 

1.0227787 

U0126725 

Lk0559915 

1.0155451 

1 • C 1577 09 
1106297CL 
1.C0C00C0 

1.0190582 

1.0713225 

1.0225396 

1.0730750 

1.0262555 

1.0694904 

1.0000000 

1.0380324 

1.0673455 

L. 0031565 
1.0432687 
1 .0540738 

1.0063769 

1.0489454 

1.0419746 

1.0096875 

1.0551352 

1.0315569 

1.0131170 

1.0619401 

1.0226102 

1.0166956 
1'. 0695 1 86 
1.0140210 

UC2C4558 
1.0701454 
l jCOOOOCO 

U 0244317 
1 k 0883560 

1.0286596 

1.0070501 

1.0331788 

1.0800652 

1.0000000 

1.0380453 

1.0660695 

1.0031588 
I .0432833 
1.0513721 

1 • 00638 06 
1.04896 16 
1.0388741 

1.009692? 

1.0551533 

1.02R5358 

1.0131226 

1.0619603 

1.0198166 

UC167C24 

1.0695411 

1.0114568 

1.C204636 
li07817C6 
1 JCOCOOCO 

1 k 0244407 
1.0883842 

1.0286698 
1L 1014070 

1.0331904 

1-0825620 

1.0000000 

1.0379359 

1.0620809 

1.0031520 

1.0431549 

1.0469301 

1.0063654 
1.0488087 
1.03464 15 

1.0096678 

1.0549655 

1.0248376 

1.0130eB3 

1.0617169 

1.0167609 

1 • C 166577 
1 k 069 1 995 
llC 092107 

l i C2C4080 
1.0776355 
UCOCOOCO 

1*0243734 

1.0874012 

1.0285900 

U0990601 

1.0330967 

1.0798632 

1.0000000 

1.0377151 

1.0569594 

1.0031364 

1.0428975 

1.0418271 

1.0063323 

1.0485048 

1.0301080 

1.0096159 

1.0545968 

1.0210987 

1.0130167 

1.0612484 

1.0138806 

1 . C 165653 
1 k 06856 1 2 
11C0734B0 

1*0202938 

1.0766840 

UOCCOOCO 

U 0242359 
1.0858363 

l .0284276 
1 k 0962779 

1.0329069 

1.0755923 

1.0000000 

1.0373934 

1.0511087 

1.0031123 

1.0425244 

1.0363431 

1.0062822 

1.0480680 

1.0254881 

1.0095383 

1.0540743 

1.0174692 

1. 0129104 
1.0605982 
1.0112392 

1.0164288 

1*0677023 

1*0057842 

110201255 
1 • 07546C 0 
l.COOOOCO 

I * 0240340 
1.0839591 

1.0281895 

1.0932962 

1.0326293 

1.0704663 

l. 0000000 
1.0369805 
1.044507P. 

1 .0030801 
1.042048? 
1 .0305263 

1.0062160 

1.0475159 

1.0208419 

1.0094366 

1.0534241 

1.0139848 

1.0127716 
1 .0598090 
1.0088239 

1.0162513 
1*066696 l 
lk 0044 4 34 

1.C199073 

UC740886 

UCCCCOCO 

lk 02 377 26 
1*0819564 

1.0278819 

1.0902662 

U0322717 
1.064404 1 

1.0000000 

1.0364847 

1.0369617 

1.0030402 
1 .04 14 BOO 
1.0243782 

1.0061346 

1.0468638 

1.0162090 

1.0093122 

1.0526691 

1.0106655 

1.0126C25 

1.0589176 

1.0066185 

1.0160356 

1.C656053 

1*0032769 

1.0196426 
UC7267 52 
UCCCCOCO 

1*0234562 

UO79970O 

1 * 0275104 
1.0872903 

1 .03 10408 
1 .0568409 

1.0000000 

1.0359137 

1.0262517 

1.0029929 
1.0408292 
1 .0179571 

1.0060388 

1.0461241 

1.0116619 

1.0091665 

1.0518272 

1.0075462 

1.0124C51 

1.0579523 

1.0046197 

1.0157844 

1.0644781 

1.0022578 

1. 0193351 
1 i 07 1 30 72 
l.COCOOCO 

1*0230091 
I. 0781744 

1.0270800 

1.0844425 

1.0313428 

1.0465251 

l. 0000000 
1.0352737 
1.0184282 

1.0029383 

1.0401032 

1.0114732 

1.0059292 

1.0453059 

1.0073285 

1.009C007 

1.0509103 

1.0046834 

1.0121813 

1.0569295 

1.0028384 

1 kO 1 55002 
lk 06 33423 
U0013743 

1.0189876 
1 4 07004 7 7 
UCOCOOCO 

1 . C226749 
1.0767494 

1.0265953 

1.0535081 

1.0307830 

1.0308121 

l. 0000000 
1-0156181 
1.0086674 

V .0012771 
1.0177804 
l .0053487 

1.00259 L7 
1.0201151 
1.0033845 

1.0039458 

1.0226350 

1-0021463 

1.0053504 

1.0253430 

1.0012926 

1.00b8186 

U0282058 

1*0006221 

l *008364 1 
1*0310146 
UCCCCOCO 

1.0100012 

1.0324798 

1.0117452 

1.0238216 

1.0136119 

1.0144440 

1.0000000 

1.0000000 

1.0000000 

1.0000000 
1.0000000 
1. 0000000 

l.OOOOOCO 

1.0000000 

l.OOOOOCO 

l-OCOOOOO 

1.0000000 

1.0000000 

l.OOOOCOO 

1.0000000 

l.CCOOCOO 

1. OCOOOOQ 

uocooooo 

i.cccoooo 

l.COOOOCO 

l.OOOOOCO 

UCCCCOCO 

UOCOOOOO 
L. OCOOOOQ 

l.COOOOOO 

i.coooooo 

l.OCCCOOO 

UOOOOOOO 


FINAL RESULTS FOR CASE 2 

FORCE/AREA = 0. 3056596E-0 1C0EF OF CENTER OF PRESSURE IN PERCENTAGE OF SIDE DIMENSIONS = I 0.4486157E 00, 0.5691036E 00). 

FLOW PER UNIT LENGTH IN X AT ENT. AND EXIT = (-0 . 205452 1 E-0 1 , 7k l 7 16304E-03 ) FLOW PER U. L. IN Y = I -1.7507014E 00) 

♦01* EXIT IN TRIPAD 
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